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PREFACE 


Exploratory modeling is presented for likely existing 
extended hydrogen atmospheres of the outer satellites 
Amalthea, Ganymede, Callisto ajid Titan. Primary emphasis 
is placed upon describing the spatial structure and Lyman-a 
intensity of these extended satellite atmospheres for a 
range of probable hydrogen atom emission and lifetime con- 
ditions. The likelyhood of detection of these extended 
atmospheres by the Voyager and Pioneer 11 spacecrafts is 
assessed. The initial spatial distribution of the ions 
created by atoms lost from the se extended hydrogen atmos- 
pheres and from the extended sodium atmosphere of lo, 
bec,^se of" ionization processes, is also presented. 

The present study provides a valuable introduction and 
theoretical description of extended satellite atmospheres, 
in addition to presenting many interesting and useful 
exploratory and prototype model calculations. Data for 
extended satellite atmospheres, obtained from both earth 
based observations and spacecraft measurements, can be 
analyzed to yield important information about the outer 
satellites and their planetary interactions. Specific 
information can be deduced about the nature of the local 
satellite atmospheres, the character of the electromagnetic 
iriteraction of satellires with their planetary magnetos- 
pheres, and the energy density and spatial distribution of 
charged particles in the magnetospheres. At present very 
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little data for extended hydrogen atmospheres has been 
accumulated. In anticipation of data to be obtained in the 
near future, exploratory modeling has been undertaken. 

Exploratory modeling results for most probable atom 
emission and lifetime conditions indicate that the extended 
hydrogen atmosphere of Amalthea should form a tightly-bound 
partial toroidal -shaped cloud about its planet. In contrast, 
Ganymede, Callisto and Titan should form rather large, 
complete and nearly symmetric toroidal-shaped clouds. The 
palatial toroidal hydrogen cloud of Amalthea, tentatively 
detected by the UV instrument aboard the Pioneer 10 space- 
craft, could be explained if the satellite atom escape flux 
were of order 10 cm sec . Model results suggest that 
the Voyager spacecraft should be able to detect the Lyman-a 
emission from the extended hydrogen clouds of Ganymede, 
Callisto and Titan, assuming the most probable atom emission 
and lifetime conditions. The Pioneer 11 spacecraft should 
also be able to detect the hydrogen cloud of Titan. Voyager 
detection of the Ganymede and Titan clouds under less favor- 
able conditions is also possible. 
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CHAPTER I 


INTRODUCTION 


Interest in understanding the atmospheres of the outer 
satellites, their extended nature and their planetary inter- 
actions, has grown significantly in recent years. Both earth 
based astronomy and spacecraft measurements have made important 
contributions. These measurements and others to be obtained 
in the near future from the Pioneer, Voyager and Galileo space- 
missions, provide a sufficiently broad data base to 
seriously constrain theoretical models and to facilitate def- 
inite physical interpretation. In this report, modeling 
focused primarily upon providing a framework for understanding 
the extended nature of outer satellite atmospheres is presented 
and applied to several satelJ.i.tes of Jupiter and Saturn. 

Toe concept of an extended satellite atmosphere is rela- 
tively new and deserves some introductory comm.ents. Briefly, 
an extended satellite atmosphere is that part of a satellite 
atmosphere which extends spatially beyond the gravitational 
control of the satellite. This situation is illustrated in 
Figure 1. Here is depicted a satellite, surrounded by a 
local or bound atmosphere, and the satellite's gravitational 
sphere of influence. Inside this sphere of influence, the 
satellite exerts more gravitational control over atmospheric 
gases than does its parent planet. If a gravitational escape 
mechanism exists, which is sufficiently energetic to allow at 
least some atmospheric gases to escape the sphere , an extended 
satellite atmosphere is created. Because of the relatively 

1 : 


Escape 



Figure 1 

Creation of an Extended Satellite Atmosphere 
See the text for discussion. 



small escape velocities of the outer satellites (between 2 
and 3 km/sec) , such gravitational escape of satellite gases 
is not unlikely. 

Gases entering the extended atmosphere in this manner 
arc basically collisionless and follow individual trajectories 
in the circumplanetary space. The spatial volume filled by 
the extended atmosphere depends upon the emerging atom veloc- 
ities and the survival lifetime of these atoms in the planetary 
environment. Cloud atoms are lost from the extended atmos- 
phere upon ionization. For short lifetimes, gases will be 
relatively confined to the near satellite environment, filling 
only a partial toroidal-shaped volume. For long lifetimes, 
however, gases will fill a complete toroidal-shaped vo3.ume, 
extending all the way around the central planet. 

Major emphasis in this report has been placed upon pro- 
viding exploratory and prototype models for extended hydrogen 
atmospheres of outer satellites which arc likely to exist and 
be discovered by the Pioneer, Voyager and Galileo outer plan- 
etary missions of NASA. Observational evidence for extended 
outer satellite atmospheres is reviewed in Chapter II. A 
theoretical framework for describing extended atmospheres is 
presented in Chapter III and the model to be used to calculate 
its detailed spatial density described. Exploratory modeling 
results for extended hydrogen atmospheres for the satellites 
Amalthea, Ganymede, Callisto and Titan are given in Chapter IV. 
Calculations describing t.ie initial ion sources created by 
these extended hydrogen atmospheres and by the sodium cloud 
of lo are also presented. Progress in modeling lo's sodium 
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line profile data 
Concluding remarks 


IS reported xn the last section of Chapter IV 
are given in Chapter V and complete the 


report . 
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CHAPTER II 


OBSERVATIONAL EVIDENCE FOR EXTENDED SATELLITE ATMOSPHERES 

Several satellites in the outer solar system - Amalthea, 
the four Galilean satellites and Titan - have, or are likely 
to have, extended atmospheres of gases, forming partial or 
complete toroidal-shaped clouds about their parent planet. 

Our present knowledge of the existance of bound and extended 
atmospheres for these six satellites is summarized in Table 1 
Discussion will be directed primarily to extended satellite 
atmospheres of atomic hydrogen gas needed for prototype and 
exploratory modeling undertaken in later chapters. 

2.1 Extended Hydrogen Atmospheres 

Extended hydrogen atmospheres, tentatively detected for 
Amalthea (JV) and definitely observed for lo (JI) with the 
ultraviolet instrumentation on the Pioneer 10 spaceeraft, 
are discussed by Judge et al. (1976) . The average brightness 
approximate vertical thickness and angular extent of these 
clouds, centered on the satellite and measured along its or- 
bital track, are given in Table 2 . Amalthea, with such a 
small escape velocity, cannot retain an atmosphere even for 
temperatures as low as ten degrees Kelvin. lo, on the other 
hand, requires a temperature of several hundred degrees in 
order for escape of hydrogen to occur . 

The other Galilean satellites Europa (JII) , Ganymede 
(Jill) and Callisto (JIV) are also likely candidates for ex- 
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Table 1 Summary of Information on the Existence of Extended Satellite Atmospheres 


Satellite 

Escape 

Velocity (km/sec ) 

Bound 

Atmosphere 

Extended 

Atmosphere 

likely: H (tentative 
detection) 

Amalthea 

0.155 

not likely 

lo 

2.56 

yes 

yes: H, Na, K 

Europe 

2.06 

probable 

yes: 0 suggested 

Ganymede 

2.75 

probable 

likely; 0, H, 

Callisto 

2.39 

probable 

likely: 0, H, H 2 

Titan 

2.73 

yes(GH^,H2) 

yes: H (tentative 
detection) 
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Table 2 Hydrogen 

Satellite 

Amal thea 
lo 


Toroidal Clouds Detected 


by the UV Photometer of Pioneer 10 


Average 

Brightness 

(RaylelghsV 

Vertical 
Thickness 
iJuplter RadH ) 

Orbital 
Ex tent 
(deg) 

100 

-1 

J 



•‘•109 

300 

- 1 



X 

- 60 
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tended hydrogen gas clouds, a source of H atoms may be 
supplxed to the cloud by water, frost or ice, detected on 
these satellites' surfaces. This might occur through dis- 
socxative surface sputtering processes by Jovian magnetos- 
pheric particles for the satellite Europa (Wu et al., 1978) , 
or through solar evaporation and photolysis for the satel- 
Ixtes Ganymede and Callisto (Yung and McElroy, 1977) , fol- 
lowed by gravitational escape. A possible detection of a 
weak hydrogen cloud for Europa may be present in the more 
recent reduction of Pioneer 10 ultraviolet photometer data 
(Wu et al., 1978) using improved procedures, although emis- 
sion from oxygen atoms may account for the complete signal 
observed. No such clouds have been detected for Ganymede 
or Callxsto in Pioneer 10 and 11 spacecraft data, although 
they may well exist, having an intensity below the sensiti- 
vi ty level of the UV instrument. 

A hydrogen toroidal cloud for Titan was first suggested 
by McDonough and Brice (1973) as a possible recycling mech- 
anism to moderate the large amount of hydrogen thought to be 
escaping the satellite. Estimates of the escape flux for 
hydrogen have' been somewhat reduced by more recent upper limits 
of the amount of present in the atmosphere of Titan 
(Mtihch et al.,1977), but a modest toroidal cloud of H atoms 
should still be expected at the satellite orbit. Tentative 
detection of this hydrogen cloud from the Copernicus satellite 
has been reported by Barker (1977) / 

E St xma ted values of the atomic hydrogen escape flux for 

each of the six satellites of Table 1 are given in Table 3 


8 



Table 3 Escape Flux of Hydrogen Atoms from the Outer Satellites 


S a t e 1 1 i t e 

Escape Flux 
(cm sec" ) 

Amalthea 

? 

lo 

O 

O 

Europa 

110® 

Ganymede 

'''5 X 10^ 

Callisto 

? 

Titan 

10^-4x10^° 


Commen t 


Reference 


possible detection, Pioneer 10 data 
calculated from Pioneer 10 data 
estimated from Pioneer 10 data 
plausible model calculation 
tnay be similar to Ganymede 
estimated from limited data 


Judge et al. (1976) 
Carlson and Judge (1974) 
Wu et al. (1978) 

Yung and HcElroy (1977) 

Tabari4 (1974) 
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I 

along with their references. An estimated escape flux for 
Amalthea is unavailable. Values for lo and Europa are inferred 
from Pioneer 10 data, whereas values for Ganymede and Titan 
are estimated from plausible model calculations. A likely 
upper limit to the hydrogen escape flux for Callisto is pro- 
vided by the Ganymede value. 

Hydrogen atoms escaping the satellites will remain part 
of the neutral extended atmosphere until ionized and then are 
swept away by the planetary magnetosphere or by the solar 
wind. For the five satellites of Jupiter under consideration, 
change exchange v/ith magnetospheric protons is thought to be 
the dominant lifetime mechanism for H atoms. Estimated life- 
time values for this process are given in Table 4 and are con- 
siderably smaller than the H atom photoionization lifetime of 

g 

about 3-6 X 10 sec. For Saturn, no measurement s are avail- 
able for magnetospheric proton fluxes at the orbit of Titan. 

If the satellite is within the planetary magnetosphere, the 
one gauss surface field of Saturn implies a proton number 
■density at Titan's orbit of about 1.7 cm"^ moving with the 
co-rotational field speed (Si scoe, 1978). This gives the 
hydrogen atom lifetime of about 8 .4 x 10^ sec shown in Table 4 
and is adopted for latter use. in the less likely event that 
Titan lies outside the Saturn magnetosphere, solar wind impact 

ionization and photoionization would result in H atom life- 
times of 2 X 10 sec and 1-2 x 10^ sec respectively. 
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Table . Estinatea biletin,e of Hydrogen in the Extended Satellite atmospheres 


Satellite 

Protcn Nuniber 
Density 

(on 

Plasma 

Tetperature 

(ev) 

Proton 

Flux 

(on ^sec~^) 

Change Exchange 
Lifetime for H 

fspr*) 

Plasma 

Conditions 

Analthea 

100 

^r. * 

100 

1.6 X 10^ 

1.8 X 10^ 

50 

Frank et al. (1976) 


:55ooo 

2 

1 X 10^° 

2.5 X 10^ 

7 

Brown (1978) 

lo 

100 

100 

1.6 X 10^ 

1.8 X 10^ 

50 

Frank et al. (1976) 


250 

13+ 

1.4 X 10^ 

2.0 X 10^ 

56 

Siscoe (1978) 

Europa 

12 

400 

3.8x10® 

7.6 X 10^ 

200 

Frank et al. (1976) 

Ganymede 

1 

400 

3.1x10^ 

9.2 X 10® 

2500 

Frank et al. (1976) 

Callisto 

1 







X 

400 

3.1 X 10 ' 

9.2 X 10® 

2500 

Frank et al. (1976) 

Titan 

1 7 







X • / 

164 

3.4 X 10 ' 

8.4 X 10® 

2300 

Sisooe (1978) 

+ assuming protons co-rotating vdth the magnetosphere 





value, calculated by BrartTi for elecfcronc? 

/ ^ is ass»ed lor pylons at Malthea’s 

orbit (2.54 Jupiter radius) as an upper limit value. 
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2.2 Sodium Cloud of lo 


In addition to th© hydrogen atoms, lo has extended atmos- 
pheres of sodium and potassium atoms, both detected from ground 
based observations. Diverse observations of the sodium cloud 
since its discovery (Brown, 1974) have produced a large and 
rapidly growing data base, sufficiently broad to seriously 
constrain theoretical models and to facilitate definite 
physical interpretation. The potassium cloud, in contrast, 
was discovered more recently (Traf ton, 1975; Trarger et al . , 
1976; Miinch et al. , 1976). To date, very little data has been 
acguired and interpretative modeling is not yet warranted. 

Observations of the sodium cloud basically are of two 
types, (1) spatial cloud data, describing the intensity of 
the sodium D- line emission in the vicinity of lo and (2) line 
profile data, giving the doppler signiture of the emitting 
cloud atoms near the satellite. Data of the first type has 
been successfully modeled by Smyth and McElroy (1978) who 
compared their model calculations with the set of observed 
two-dimensional images reported by Murcray and Goody (1978). 
These results suggest that sodiuiri atoms originate from the 

inner ■’.emi sphere of lo with an average flux of about 10^ 

. -2 -1 . , 
atoms cm sec and an initial mean velocity of 2.6 km sec~"^. 

Atoms appear to be removed from the cloud by electron impact 

ionization with an effective lifetime of between 15 and 20 

hours. Data of the second type, the line profile observations, 

has received less attention. Preliminary modeling efforts 

(Carlson et al . , 1978 ; Smyth and McElroy, 1977) need to be 
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continued and a comprehenaive program undertaken. Our develop- 
ment of such a modeling program is already in progress and 
will be discussed in a later chapter. 
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CHAPTER III 


MODELS FOR EXTENDED SATELLITE ATMOSPHERES 

A comprehensive physical description for extended satel- 
lite atmospheres is presented in this chapter. A simple 
theory, useful in initially estimating the spatial size of 
extended satellite gas clouds, is summarized and contrasted 
with a more exact physical description used in our cloud model. 
The simple theory is applied to the one existing and five 
likely existing hydrogen atmospheres discussed in the previous 
chapter and overall spatial dimensions of these clouds are 
estimated. A brief description of the cloud model, capable 
of calculating the detailed intensity and density variations 
or line profile shapes of extended satellite atmospheres, is 
Lhon presented. 

3.1 Physical Description of Extended Satellite Atmospheres 

The atmosphere of an outer satellite may be divided into 
two parts: a portion that is gravitationally bound to the 

satellite and a portion that is gravitationally unbound, the 
so-called extended satellite atmosphere. The spatial boundary 
between these two portions of the atmosphere occurs where the 
gravitational force fields of the planet and satellite on an 
atom are comparable. The surface is effectively a sphere 
centered on the satellite with a radius r equal to the 

Xj 

distance between the satellite and the near collinear Lagrange 
point. An approximate expression for the Lagrange radius is 
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^ 


3 (1-y ) 


1/3 


(3.1.1) 


where a is the orbital radius of the satellite and „ is the 

ratio of the satellite mass to the sum of the planet and 

satellite masses Values of r^ as well as several other orbital 

and physical properties of the six outer satellites are given 
in Table 5 . 

an atom emitted from the surface or exosphere of a satel- 
lite will be confined to the Lagrange sphere unless it has a 
velocity nearly egual to the satellite escape velocity (see 
for values). The situation, including the effect of 
duplter gravitational field, is illustrated in Figure 2 for 
Ganymede, where an atom orbit, initially emitted from the 
exosphere, is shown as a function of increasing initial 
velocity. The escape velocity is reduced from 2.6 km/sec, 
the two-body value, to about 2.45 km/sec because of the 
gravxtational attraction of Jupiter. 


3.1.1 Angular Extent of the Toroidal Clouds 

The orbital path of atoms escaping the Lagrange sphere 
will be controlled largely by the gravity of the planet. 

Those escaping atom orbits with Initial emission velocities 
less than about (/J -i), „here v^ is the orbital speed of 
the satellite (see Table 5), will have elliptical orbits about 
the planet. For isotropic emission of many such orbits, atoms 
will diffuse both ahead and behi.nd the satellite. An approx- 
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Table 5 Physical and Orbital Properties of the Outer Satellites 


Satellite 

Satellite 

Radius 

(km) 

Orbital 
Radius 
( 1Q5 km) 

Orbital 

Speed 

(km/sec) 

Orbital 

Period 

(hr) 

Satellite 
Mass^ 
(1022 kg) 

Lagrange 

(km) 

Radius r 
(Satellite 
radii) 

Exosphere 

Radlus^^ 

(km) 

Escape 

Velocity* 

(km/sec) 

Amalthea 

120 

1.813 

26.47 

11.96 

2.17 X 10~^ 

283 

2.36 

120 

0.155 

lo 

1820 

4.216 

17.33 

42.46 

8.91 

10,542 

5.79 

2600 

2.14 

Europa 

1525 

6.709 

13.74 

85.23 

4.87 

13,717 

9.00 

2000 

1.80 

GanjTneda 

2635 

10.700 

10.88 

171.71 

14.90 

31,759 

12.05 

2935 

2.60 

Gal 11s to 

2500 

18.800 

8.19 

400.54 

10.74 

50,031 

20.01 

2800 

2.26 

Titan 

2500 

12.220 

5.57 

382. 6 «J 

14.01 

53,106 

21.24 

4-10 x 10^ 

2.16-1.37 


Mass of Amalthea calculated assuming a density of 3 gm cm and a radius of 120 km. 
Estimated values assuming a bound atmosphere except for Amalthea. 

* Two-body escape velocities computed for the exosphere radius Indicated. 
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BALLISTIC ORBtTS OF GANYMEDE (Gonymede's Fiona) 


TOWARD JUPITER 


t 




VELOCITY MAGNITUDE (km/stc) 2.44 2.4425 


Figure 2 




Ballistic and Escape Orbits of Atons Qnitted frcm Ganiroede. Orbits of an atom 
initially ardtted frori Ganynede's exo^here, and directed radially av.ay frcsn 
Jupiter in the orbital plane of the satellite, are shown for increasing values 
of the atari initial velocity magnitude. The orbits are shown as seen by an 
observer on Ganymede and the tick marks occur at 2 hr. intervals along the orbit. 
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imate value for the forward and backward angular diffusion 
velocity of the resulting cloud relative to the satellite, 

estimated by using the expressions of Fang 
et al. (1976), which assumes a massless point satellite with 
an isotropic Maxwell Boltzmann initial speed distribution. 
Results of this calculation for the six outer satellites 
of interest are given in Figure 3 through Figure 8 and compared 
with the approximate expression for the angular diffusion 
velocity 


a = 1 1/2 

a irm 


(3.1.2) 


where T is the temperature and m is the mass of the emitted 
atom. The emission temperaturo T roughly oorresponds to an 
atom emerging from the Lagrange sphere with an average velo-^ 
city V, relative to the satellite, given by 

V = (^) = 7.28 X 10"^ km/sec , (3.1.3) 


where T is in °K and A is the atom msss in AMU. 

The angular extent of the forward cloud 0„ and backward 
cloud 0g about the planet, relative to the satellite location, 
is proportional to the lifetime t of the neutral cloud atoms 
in the planetary environment 


®F ,B “ f , B 


(3.1.4) 
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Figure 3 


--^o- -uplte. o. 

fip and Q are both compared with the simnlA^^^^ velocitie 

temperatgre to atomic weight ratio ^ expression with the 

dependent variable. emitted atoms as the 
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Figure 4 

Forward and Backward Angular Diffusion Velocities 
about Jupxter of the Gas Cloud Emitted by lo 
Same description as Fig. 3. ^ 


AVERAGE ANGULAR DISPLACEMENT FROM lo PER lo PERIOD (deg) 



AVERAGE ANGULAR DIFFUSION VELOCITY (deg/hr) 



AVERAGE ANGULAR DISPLACEMENT FROM EUROPA PER EUROPA PERIOD (deg) 




SATELLITE : GANYMPDF 



temperature/atomic weight {®K/AMU) 


Figure 6 


Forward and Backward 
about Jupiter of the 
Same description as 


Angular Diffusion Velocities 

Gas Clouds Emitted by Ganymede. 
Fxg. 3. 


AVERAGE ANGULAR DIFFUSION VELOCITY (deg/hr) 



— f Backward Angular Diffusion Velocities 
the Gas Clouds Emitted by Cal- 
iisto. Same description as Fig. 3. 



AVERAGE ANGULAR DIFFUSION VELOCITY (deg/hf) 



Figure 8 

Forward and Backward Angular Diffusion Velocities 
about Saturn of the Gas Clouds Emitted by Titan. 
Same description as Fig. 3. 


AVERAGE ANGULAR OISPLACEHENT FROM TITAN PER TITAN PERIOD (deg) 



The forward and backward cloud will diffuse through an angular 
distance of 360° with respect to the satellite in a time inter- 
val Tp of order 


with a cylindrical symmetrical toroidal cloud forming if the 

lifetime x is somewhat in excess of x„^ . Values of x ^ 

f B F , B 

for atoms with an emission temperature to mass ratio of 
200° K/AMU are given in Table 6 for the six satellites and 
corresponds to hydrogen atoms escaping from the satellite 
Lagrange sphere with a mean velocity of 1.03 km/sec. 

The gravitational field of the satellite naturally com- 
plicates this simple description near the Lagrange sphere. 
This is illustrated for the satellite Ganymede for both 
molecular and atomic hydrogen toroidal clouds. Thermal 
emission is assumed and the time evolution of the cloud is 
examined by solving the three-body gravitational problem. 

Using the Ganymede atomspheric model of Yung and McElroy 
(1977), an exospheric radius of 2935 km and a temperature of 
14 0 K are adopted, where II and II 2 have an assumed Maxwell 
Boltzmann speed distribution ^ Restricting attention to 
particle orbits in the satellite plane only, the time evolu- 
tion of the column density of the two toroidal clouds is 
shown in Figure 9 and Figure lO. Ganymede is positioned at 
the elongation point of its orbit in the Figures. A symmetric 
toroidal cloud is established for atom flight times of order 
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Table 


6 Angular Diffusion Times 

F 

and Tg for Outer Satellite Clouds'^ 

Satellite 


D 

(hr) 

D 

ihrl 

Amalthea 


106 

101 

lo 


248 

233 

Europa 


402 

371 

Gan>Tnede 


654 

593 

Calllsto 


1189 

1074 

Titan 


835 

808 

+ Calculated assuming 

an atom emmlssion temperature of -200® 

K/AMU 
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Figure 9 

Time Development of the H Atom Toroidal^ Cloud in Ganvmede ' s 
orbital Plane. Each columji density profile is calculated as 
described in the text. Ganymede's location at the elonaation 
point of the orbit is indicated . 
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COLUMN DENSITY (ARBITRARY UNITS) 



DISPLACEMENT FROM GANYMEDE (6xlO‘*km) 


C^MEDE'sqRBIT 
30 40 50 


Figure 10 

Time Development of the Molecular Toroidal Cloud 
in Ganymede ' s Orbital Plane. Same descriotion as 



I 


of 800 - 1200 hours, excluding the immediate vicinity of 
Ganymede. The column density near Ganymede is larger because 
84% of the atoms and 98% of the H 2 molecules emitted per unit 
txme are gravitationally confined to the Lagrange sphere. 

3.1.2 Vertical Extent of the Toroidal Clouds 

The height of the toroidal cloud above the satellite 
plane H xs only a function of the initial emission velocity 
xf the lifetime is longer than about one-fourth of the satel- 
Ixte perxod, a time required for an atom orbit to reach max- 
imum vertical excursion. A simple expression for the height 
of the cloud directly above the satellite orbit is given by 


where the latter expression follows from the definition (3.1.3) 
3.1.3 Radial Extent of the Toroidal Clouds 

For a cordinate frame centered on the planet, the inner 
radial boundary r_ and the outer radial boundary r^ of the 
gas toroidal cloud are given by the simple two-body formula 

''t = — — V ^ — • (3.1.7) 

2 ( — ^ ) ^ -1 


= 3- -V2 (3.1.6) 

■ Vg \m ’ 
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This is valid if the atom lifetime t is equal to or larger 
than half the atom orbit period T_j_, given by 



where is the period of the satellite. 

3.2 Estimated Spatial Extent of the Extended Hydrogen Atmospheres 

In order that a satellite have an extended hydrogen atmos- 
phere , it must have 

(1) a source of hydrogen atoms 

(2) a sufficiently energetic atmospheric process for 
gravitational escape of these atoms, and 

(3) a planetary environment providing a sufficiently 
long lifetime for escaping atoms so they can popu- 
late a spatially extended volume. 

The shape and size of the gas clouds depends critically 
upon the emission conditions (2) and the lifetime (3), 
whereas the density within this volume depends upon the mag- 
nitude of the escaping atom flux. Estimated values of the 
hydrogen escape flux for the outer satellites are summarized 
in Table 3 . These flux values suggest that detection of 
extended hydrogen toroidal clouds for Europa, Ganymede and 
possibly Callisto (if similar to Ganymede) should be of 
order 100 times more difficult than for lo. Titan's cloud, 
on the other hand, should be more comparable to lo's hydro- 
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gen cloud and may well be more easily detected. 

In What follows, it will be assumed that the first con- 
dition above is satisfied. Lifetime of hydrogen atoms in the 
planetary environment, needed in evaluating the third con- 
dition, are given in Table 4. what then remains, in order 
to estimate the spatial extent of the toroidal hydrogen atmos- 
phere, is specification of the emission conditions for H atoms. 

3.2.1 Emission Conditions 


The emission conditions for H atoms - the speed, disper- 
sion and angular character of the initial atom velocity vec- 
tors over the exosphere surface of the satellite - are not 
known for the satellites, m the absence of such information, 
attention will be directed to isotropic emission assumincj a 
thermal-like escape mechanism. Thermal -like escape here will 
describe atoms which have an emission velocity distribution 
characterized by an exospheric temperature. For outer satel- 
lites, a value 250" K corresponds to hydrogen atoms with a 
most protable exospheric emission velocity of 2.0 km/ sec, 
assuming the simplest case of a Maxwellian distribution. 

Thrs emission velocity is comparable to the escape velo- 
cities of the outer satellites (see Table 5) . The bulk of 
such H atoms, emitted from the exosphere and emerging from 
the Lagrange sphere, will therefore have an Initial speed 
restricted to the tail of the thermal distribution and an 

speed of order 1 km/sec. This exit velocity corresponds 
for hydrogen to a two-body atom emission temperature, used 
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earlier in Figure 3 through Figure 8, of between 100 and 
200° K/AMU. 

3.2.2 Spatial Extent Estimates ' • 


The angular extent of the extended hydrogen oloude of 
the Six outer satellites can now be estimated fro„ the ex- 
pression (3.1.4) using the angular diffusion velocities of 
Prgure 3 through Figure 8 and the estimated lifetime values 
of Table 4. Results are given in Table 7 for atom emission 
temperature of 100, 200 and 500- K/h«o, and correspond to a 
velocity V for hydrogen atoms emerging from the bagrange 

spheres of about 0.73, 1.03 and 1 Vm/ 

and 1.63 km/sec. The lower two 

velocities represent thermal emission as discussed earlier 
The highest velocity case is more appropriate for „on-thor- 
mal emission such as surface sputtering or electromagnetic 
Oriven escape processes. Observations for lo and Amalthca 
(Table 2) suggest that thermal emission may be more appro- 
priate, resulting in partial toroidal clouds for hmalthea, 

TO and Kuropa, and somewhat well developed complete toroidal 

Clouds for Ganymede, Callisto and Titan. 

The vertical and radial extent of these clouds may be 
calculated using the expressions (3.1.6) and ( 3.1 7 ) 

Kesults are given in Table 8 for an atom emission tempera- 
ture of 200. K/A,,o. The estimated total height of 0.7 
Jupiter radii for the lo clonr^ to • 

Cloud IS in agreeinent With the Ob- 

servations of Table 2. Note that in Table 8 the individual 
clouds for the Galilean satellites overlap radially, so 
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Table 7 Estimated Forward and Backward Angular Extent of the Hydrogen Toroidal Clouds 



100® 

K/AMU 

200® 

K/AMU 

500® 

K/AMU 



0 

0 

0 

0 

0 

0 

Hydrogen 


F 

B 

F 

B 

F 

B 

Lifetime 

Satellite 

(deg) 

(deg) 

(deg) 

(deg) 

(deg) 

(deg) 

(hr) 

Amalthea 

120 

130 

170 

180 

270 

280 

50 


17 

18 

24 

25 

37 

40 

7 

lo 

50 

50 

70 

80 

110 

120 

50 


60 

60 

80 

90 

120 

140 

56 

Europa 

130 

140 

180 

190 

270 

300 

200 

Ganymede 

1000 

1100 

1400 

1500 

2000 

2200 

2500 

Callisto 

560 

600 

760 

840 

1100 

1100 

2500 

Titan 

750 

830 

990 

1030 

1340 

1180 

2300 


Table 8 Estimated Vertical and Radial Extent of the Hydrogen Toroidal Clouds + 


Satellite 



Satellite 


H 

r 



(Planet radii)* 

(Planet radii) 

Location 
(Planet radii) 

^+ 

(Planet radii) 



Amalthea 

+0.099 

2.18 

2.54 

2.98 

lo 

+0.351 

4.69 

5.91 

7.56 

Europa 

+0.704 

7.03 

9.40 

12.86 

Ganymede 

+1.003 

10.42 

15.00 

22.41 

Callisto 

+3.309 

16.31 

26.35 

45.54 

Titan 

+3.74 

10.08 

20.23 

47.58 


+ Calculated assuming an- atom emission temperature of 200*^ K/AMU 
* Planet Radius, 7.135 x lo'* for Jupiter and 6.04 x 10^ km for Saturn 


pro- 


that if each one exists, a spatially extended source of 

tons for the Jovian magnetosphere Is provided in the radial 
annulus from 4.7 to 45.5. 

3 • 3 Cloud Model 

In the cloud model, atoms are released radially from a 
satellxte's exobase with some specified distributions of 
speeds and emission flux. Atom trajectories are modeled by 
numerically solving the restricted three -body problem as 
described by Smyth and McElroy (1977) . The contribution of 
individual atoms to the cloud density in any given volume of 
space is proportional to the initial source strength and the 
residence time appropriate to the volume element. Atoms are 
assumed emitted continuously and at a constant rate, with 
individual trajectories terminated according to some prc- 
specified lifetime, GalGUlations were performed on a CDC 
7600 computer with a typical model for the cloud assuming 

that the exobase is divided into 1298 source elements. 

For extended hydrogen atmospheres, cloud intensity 
contours are calculated assuming solar resonance scattering 
of Lynan-u radiation by the hydrogen atoms. m the case of 
sodium, the radiation intensity in the D lines is calculated 
for the optically thin limit, also assuming resonance scat- 
tering of sunlight. Proper account is taken of Doppler shifts 
introduced due to the motion of sodium atoms with respect to 
the sun and the Fraunhofer absorption feature present in the 

solar spectrum. 
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CHAPTER IV 


MODELING RESULTS 

4.1 Introduction 

Exploratory prototype modeling of hydrogen toroidal gas 

clouds for the outer satellites Amalthea, Ganymede, Callisto 

and Titan has been undertaken using our cloud model. Model 

calculations are performed for the probable atom emission 

and lifetime conditions outlined in earlier chapters. Results 

suggest than the tentative detection of a partial toroidal 

hydrogen cloud for Amalthea by the UV instrument of Pioneer 10 

*^ould be explained if the satellite atom escape flux is of 
11 —2 —1 

order 10 cm sec . Model calculations confirm the earlier 
estimated results of Chapter III that the extended hydrogen 
atmospheres of Ganymede, Cal.listo and Titan form complete 
toroidal- shaped clouds about their parent planets. The 
spatial structure, dynamics and detectability of these hydrogon 
clouds are discussed in the following section. 

The initial spatial distribution of protons created by 
the modeled extended hydrogen atmospheres of Amalthea, 

Ganymede, Callisto and Titan is calculated and presented 
in the third section. Also included is the initial spatial 
distribution of heavy ions created by the sodium cloud of lo. 
r‘inal].y, the current status of modeling the line profiles of 
sodium atoms emitting D1 and d 2 radiation from lo ' s extended 
sodium atmosphere is discussed in the last section of this 
chapter . 
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4.2 Exploratory Modeling of Extended Hydrogen Atmospheres 


The model computations reported here were performed using 
our cloud model described in section 3,3, which includes both 
the gravitational effects of the satellite and planet upon 
the extended-atmosphere atoms. For the four satellites of 
interest, the most probable values for atom emission condi- 
tions and lifetime, estimated earlier, were assumed to define 
modeling parameters. For each satellite, additional values 
for these parameters were also considered to illustrate the 
variety of results that might be expected due to uncertain- 
ties in prescribing the atom emission and lifetime conditions 
Isotropic emission of hydrogen atoms from the satellite exos- 
phere was assumed in all model calculations. 

The resulting spatial structure of each toroidal-shaped 
extended atmosphere is presented in the form of several two- 
dimensional contour plots. For the case of intensity plots, 
the contours were calculated assuming resonance scattearing 
of solar Lyraan-a radiation by the hydrogen atoms. 

4.2.1 Results for Amalthea 

Primary emphasis in modeling the extended hydrogen atmos- 
phere of Amalthea has been to estimate the atom emission 
parameters necessary to form an H cloud, with Lyman-a inten- 
sity comparable to that tentatively detected by the UV 
instrument aboard the Pioneer 10 spacecraft (Judge et al . , 
1976) . The most favorable condition for the UV instrument to 
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on spacecraft time 


detect the H cloud occurred at 18.0 hours 

day 334 GMT (see Judge et al., 1976; P-1075) when Amalthea 
was in the field of view. 

A model computation, assuming isotropic emission of H 
atoms from the satellite surface with a radial velocity of 
1 km/sec, predicts a partial toroidal -shaped hydrogen cloud 
of approximately the correct angular extent for a cloud 
lifetime of 50 hours. The column density contour plot of 
thxs partial toroidal cloud viewed from above the satellite 
plane is shown in Figure 11. This model computation further 
indicates that the signal measured at 18.0 hours on day 334 
could be produced by solar resonance scattering of Amalthea ■ s 
extended hydrogen atmosphere if the escape flux of hydrogen 
from the satellite surface were of order 10^^ atoms cm’^sec"^ 
The Lyman-a intensity distribution of this cloud for a 
10 - cm sec ^ flux is shown in Figure 12 from a Pionoor 10- 

Ixke viewing perspective nearly parallel to the satellite 
plane. The one degree width by ten degrees high aperature 
of UV instrument, assumed centered on the satellite for this 
calculation, sees all the lyman-a radiation emitted from the 
cloud Within a slice 0.52 Jupiter radii left and right of 
Amalthea. A flux of lO^^ cm"2sec-l, although large and per- 
haps improbable, is not impossible. if supplied from the 
satellite surface, for example from H 2 O ice , it would be 
equivalent to depleting only a few kilometers of the surface 
in 4 X 10^ years. if supplied by satellite neutralization 
of energetic magnetospheric protons , a large flux of protons , 
xn excess of 10 ^^ cm ^sec“^, would be required. 
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AM ALTHEA 50 hr. 



velocxty and a 50 hour u'uuex lor a i.o kin/ sec emxssion 

satellite plane. The outer ISnioS/valSrifaSnrJ/foTt'S''" "''• 



AM ALTHEA 50 hr. 



Lyitian-a Radiation from Analthea ' s Extended Hydrogen Atmosphere. 
Intensity contours / calcuXated assuming resonance scattering of 
sunlight and the model parameters indicated, are shown for the 
satellite orbit plane tilted by 3.254 degrees. 




In the event that a higher emission velocity or a longer 

lifetime is appropriate for the ertenaed hydrogen atmosphere 

of Amalthea, the resulting toroidal cloud will exist all the 

way around Jupiter. This situation is illustrated in Figure 13 

and Figure 14 where model results are shown for an emission 

velocity of 1.5 km/sec and a lifetime of 50 hours. The 

toroidal cloud, although complete is hv nr. n, 

Mprere, xs by no means symmetrical 

about Jupiter. 

4.2.2 Results from Ganymede 


Unlike Amalthea, the estimated lifetime for Ganymede -s 
hydrogen cloud atoms is sufficiently long ( 2500 hr) so as 
to ensure formation of a complete, nearly symmetric, toroidal- 
shaped gas cloud. The time evolution of this cloud is 
illustrated in Figure 15 through Figure 19, where the circular 
orbat of the satellite is viewed from above the orbit plane 
■the colu,„n density contours are calculated assuming radial 
emission of H atoms from the satellite exosphere (2935 km 

radius) with an initial velocitv ot 9-71/ 

vGxocxty of 2.7 km/see. Such atoms 

emerge from the Lagrange sphere with a velocity of about 

1 km/sec and the cloud dynamics correspond approximately to 

the 200“ K/AMH atom emission temperature results discussed 

n Chapter III. The forward and backward clouds diffuse 

enind the satellite (see Figure 15-17) with their 

outer contours (having, about 1/20 of the peak column density * 

value) meeting on the far side of the circular orbit in about 

600 hours Of flight time. For longer times (Figure 18-19) 
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AM ALTHEA 50 hr. 


-92 Royleighs 


Max Signal = 738 Rayleighs 


Model Parameters 

Emission Velocity = 1.5 km/sec 
Atom Lifetime = 50hrs. 
Hydrogen Flux =10^^ cm'^sec“’ 


Figure 14 

Lyman-a Radiation from Amalthea ' s Extended Hydrogen 
Atmosphere. Same description as Fig. 12. 



GANYMEDE 200 hr. 



Figure 15 

Cloud Dynamics for the Extended Hydrogen Atmosphere of Ganymede. Columr 

' calculated using the cloud model for an initial soeed 
and the Ixfetime indicated above, are shoxvn as viewed" fron 
above the satellite plane. The location of Ganymede, its orbital track 
and Jupiter are indicated. 



GANYMEDE 400 hr 



Figure 16 


Cloud Dynamic 
of Ganymede. 


c i extended Hydrogen ; 
Same description as Fig? 15, 



Cloud Dynamics for 
of Ganymede . Same 


the Extended Hydrogen Atmosphere 
description as Fig. 15, 
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GANYMEDE 1200 hr. 



Figure 18 

Cloud Dynamics for the Extended Hydrogen Atmosphere 
of Ganymede. Same description as Fig. 15. 
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GANYMEDE 2500 hr. 



Figure 19 

Cloud Dynamics for the Extended Hydrogen Atmosphere 
of Ganymede. Same description as Fig. 15. 
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the two clouds coalesce and produce a complete toroidal 
cloud which becomes increasingly more cylindrically symmetric. 

Voyager spacecrafts approaching Jupiter in 1979 will 
view the cloud nearly in the plane of the satellite so that 
the circular orbit appears as a thin ellipse. Results given 
in Figure 19 are also presented in Figure 20 as they would 
appear to the Voyager spacecrafts. Intensities resulting 
from resonant scattering of solar Lyman-a radiation by 
cloud atoms are indicated assuming an escape flux of 5 x 10^ 

H atoms cm ^sec ^ from the satellite. The UV instrument 
aboard the Voyager spacecrafts is capable of detecting a 5 
Rayleigh signal above the background ( 150 Rayleights) for 
a two hour period of observation (Broadfoot 1978). This 
instrument , with a rectangular viewing aperature of 0 . 1 by 
1.0 dogreos/ sees a field of view of 0 .1 by 1 . 0 diametors 
of Jupiter when the spacecraft is about 60 planetary radii 
from the satellite. The calculatod cloud intensity should 
therefore be easily detected by the UV instrument when ob- 
serving near the end of the satellite orbit. Detection be- 
comes more marginal when viewing spatial regions closer to 
Jupiter. Absence of a detectable signal will provide an 
upper limit to the hydrogen escape flux. 

If the lifetime of hydrogen atoms in the vicinity of 
Ganymede ’ s orbit were much shorter than the 2500 hour 
estimate, detection of the cloud by Voyager would be more 
difficult. This situation is illustrated in Figure 21 and 
Figure 22 where model results are shown for a lifetime of 
200 and 900 hours respectively. The emission velocity and 
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GANYMEDE 2500 hr. 



I 

Max Signal s 42.2 Rayleighs 


I 


Max Signal = 32.6 Rayleighs 


Model Parameters 

Emission Velocity = 2.7 km/sec 
Atom Lifetime = 2500 hrs. 
Hydrogen Flux = SxIO^cm'^sec’ 


rxgure 


Lyman a Radiation from Ganymede ' s Extended Hydrogen Atmosphere, 
ntensity contours, shovm lor the satellite orbit plane tilted 
by 3.2d4 degrees , v/ere calculated assuming resonance scatterina 
of sunlight and the model parameters indicated. ttering 



GANYMEDE 200 hr. 


6-3 y- 4.9 Rayleighs 



3.5 ^0.66 


t 

Max. Signal = 15 Rayleighs 

Model Paronfieters 

Emission Velocity = 2.7 km /sec 
Atom Lifetime = 200 hrs. 
Hydrogen Flux =5x10^ cm"^sec"* 


Figure 21 


Lyman-c: Radiation for Ganymede's Extended 
Atmosphere. Same description as Fig. 20. 


Hydrogen 
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GANYMEDE 900 hr. 



Max. Signal = 19 Royleiahs » . ■ 

^ Max Signal = 10 Rayleighs 


Model Parameters 

Emission Velocity = 2.7 km/sec 
Atom Lifetime ' 900 hrs. 
Hydrogen Flux = 5x10^ cm"^sec‘* 


Figure 22 


Radiation for Ganymede ' s Extended Hydrogen 
Atmosphere. Same description as Pig. 20. ^ 
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atom flux are assumed to be 2.7 km/sec and 5 x lo"^ cm'^sec"^ 
repsectively, as before. For the smaller lifetime, the cloud 
is too faint to be seen by Voyager except very near the satel- 
lite. Detection for the 900 hour lifetime case would be ex- 
pected near Gaqymede and at the two elongation points of the 
satellite orbit. 

4.2.3 Results for Callisto 

Similar to Ganymede, the estimated lifetime for Callisto 's 
hydrogen cloud atoms is sufficiently long to allow a complete 
nearly symmetric toroidal-shaped gas cloud to develop. For 
model computations, the initial velocity of H atoms emitted 
radxally from the satellite exosphere is assumed to be 2.4 
km/sec, so that atoms einercjing from the Lagrange sphere have 
a velocity Of about 1 km/sec. The flux of H atoms emitted 
by the satellite is chosen to be 5 x 10^ cm“^sec~^, the 
value suggested for Ganymede by Yung and McElroy (1977). 

This value is likely on upper bound since Callisto has less 
H 2 O ice or frost surface cover than Ganymede , and would then 
be less capable of supplying H atoms through postulated 
evaporation, photodissociation and escape processes. 

Model computations, adopting the above conditions and a 
hydrogen lifetime of 2600 hours, are shown in Figure 23 and 

illustrate the nearly symmetric structure of the cloud as seen 
from the viewing perspective of a Voyager spacecraft in Figure 24 
Here cloud intensity contours are less than 5 Rayleighs 
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CALL I STO 2600 hr. 


Viewing Aperature of Voyager 
UV Instrument (at 60 Rj) 




Figure 24 

Lyman-a Radiation for Callisto's Extended Hydrogen Atmosphere. Intensity 
contours, calculated assuming resonance scattering of sunlight and the 
indicated model parameters, are shovrn for the satellite orbit plane 
tilted by 3.254 degrees. 



(Voyager's UV instrument detection limit) everywhere except 
very near the satellite and near the two orbit elongation 
points. Intensities near the satellite are no larger than 
about 10 Rayleighs and signals at the elongation points are 
just at the 5 Rayleigh limit. In the best of circumstances 
then, detection of the extended hydrogen atmosphere of 
Callisto is only marginal in the immediate vicinity (one 
Jovian radius) of the satellite - 

The probability of detecting a hydrogen toroidal cloud 
for Callisto would be considerably improved if the lifetime 
of H atoms were actually larger than -he value adopted here. 
This would be the case if the proton flux continues to dimin- 
ish beyond Ganymede' s orbit. The lifetime adopted here 
assumes it it; constant. No data is presently available for 
the radxal var ration of thu low-enorgy proton flux in this 
spatial region. Observations of a hydrogen gas cloud for 
Callisto, significantly brighter and more symmetric than 
that shown in Figure 24 , would he consistent with such a 
radial behavior . 

4.2.4 Results for Titan 

The dynamic development of the extended hydrogen atmos- 
phere of Titan is immediate to that of Ganymede and Callisto, 
for similar emission conditions. Model calculations were 
performed assuming H atoms to have an initial velocity of 
2 . 0 km/sec at the satellite exosphere (assumed to have a 
nominal radius of 5000 km) . Atoms then emerge from the 
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Lagrange sphere with a velocity of about 1 km/ sec. Estima- 
tions of the H atom escape flux from Titan (defined with 
respect to the satellite surface) range from about 1 x 10^ 

to 4 X 10^^ cm ^sec"^ (Tabarie, 1974). A value of 1 x 10^° 
-2 -1 

cm sec has been adopted. Assuming these emission condi- 
tions and a hydrogen lifetime of 2300 hours, a column density 
contour plot of the toroidal hydrogen cloud of Titan is shown 
in Figure 25 as seen from above the satellite plane. The 
cloud has a fair amount of circular symmetry except near the 
satellite where the column density is peaked. 

The Pioneer 11 and Voyager spacecrafts will see this 
cloud more nearly in the satellite plane as shown in Figure 26 
where intensity contours, assuming solar resonance scatter- 
ing, are plotted. The cloud intensity is several hundred 
Rayleighs near the elongation points of the orbit, whore the 
geometry of the toroidal -shaped atmosphere enhances the gas 
column length and even 50 to 100 Rayleighs near Jupiter. 

Such a signal should be detected by the Pioneer 11 spacecraft 
when viewed through the 1 x 10 degree aperature of its UV 
instrument. Through the smaller aperature ( 0.1 xl degree) 
of the UV instrument on the Voyager spacecrafts, the cloud 
should be detected even if the H atom escape flux from Titan 
were reduced by a factor of 10 or 100. 

In the event that the lifetime of hydrogen atoms near 
the orbit of Titan is actually significantly smaller than 
the 2300 hour estimate, a partial toroidal cloud could exist. 

This is illustrated in Figure 27 for a lifetime of 400 hours. 
For more energetic emission conditions, the cloud could also 


Figure 25 


Model for the Extended Hydrogen Atmosphere of Titan. Column 
density contours, calculated from the cloud model for a 2.0 
km/sec emission velocity arid a 2300 hour lifetime, are shov/n 
as viewed from above the satellite plane. 



TITAN 2300 hr. 



Max. Signal - 901 Rayleighs 


Max. Signal =738 Rayleighs 


Model Parameters 

Emission Velocity = 2.0 km/sec 
Atom Lifetime = 2300 hrs. 
Hydrogen Flux =IO'Ocm'2sec'’ 


Figure 26 
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TITAN 400 hr. 



Figure 27 


Lyman-a Radiation 
Atmosphere, Same 


from Titan * s Extended Hydrogen 
description as Fig. 26. 



6 0 




be complete for this shorter lifetime, although not vertically 
and cylindrically symmetric. This case is illustrated in 
^io^ure 28 for an emission velocity of 2.6 km/sec and a life- 
time of 450 hours. 

4.3 Ion Sources Created by Extended Satellite Atmospheres 

The presence of extended satellite atmospheres of neutral 
have significant effects upon an existing planetary 
magnetosphere. Atoms lost from the toroidal atmosphere 
through ionization processes, could provide a magnetospheric 
ion source which may locally dominate other ion sources 
supplied by the planet or by the solar wind. As a first step 
in understanding this, the initial spatial distribution of 
ions created by extended satellite atmospheres will be esti- 
mated using our cloud model. Results will bo presented for 
the prototype hydrogen clouds of Amalthea, Ganymede, 

Callisto and Titan discussed in sjoction 4.2 and the sodiujii 
cloud of Jo modeled by Smyth and McElroy (1978). 

4.3.1 Proton Sources 


An average radial distribution of hydrogen atoms in an 
extended satellite cloud may be calculated by integrating 
over the angular. and vertical dimensions of the cloud. 
Dividing this distribution by the lifetime of the cloud atoms 

provides a radial distribution for the loss rate of hydrogen 
atoms. This loss rate distribution for Amalthea, Ganymede 
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TITAN 450 hr 


r35.1 




Max Signal = 234 Royleighs Max. Signal --43.3 Rayleighs 


Model Parameters 

Emission Velocity = 2.6 km /sec 
Atom Lifetime = 450 hrs. 
Hydrogen Flux =10^®cm“2sec‘^ 


Figure 28 


Lyman-a Radiation from Titan's Extended Hydrogen 
Atmosphere. Same description as Pig, 26. 
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and Callisto is shown in Pivjure 29 and corresponds to the 
modeled extended atmosphere presented, respectively in 
Figure 12, Figure 20 and Figure 24 . The loss x"ate distribu- 
tion for Titan is given in Figure 30 and corresponds to the 
model results shown in Figure 26. 

In Figure 29, the maximum hydrogen loss rate per 0.1 Jupiter 

radii for Amalthea is about a factor of 50 larger than that for 

Ganymede and about a factor of 100 larger than that for Callisto. 

For comparison, the maximum value provided by lo's hydrogen 

cloud (located at 5.9 on the horizontal scale in Figure 29) 

26 

might be of order 1 x 10 per 0.1 Jupiter radii per second, 
but detailed modeling not performed here is needed to more 
accurately specify this number. Although the loss rate per 
0.1 Jupiter radii for Anialthea and lo would be comparable 
under Lhese circumstances, the total number of hydrogen atoms 
lost per second by lo would be of order 10 times greater than 
for Amaltheci. in Figure 30, the maximum hydrogen loss rate 
per 0.1 Saturn radii per second by Titan is about 1.5 x 10^^, 
with the complete cloud losing about 7.8 x 10^^ hydrogen 
atoms per second. 

Assuming a steady cloud state, the hydrogen atom loss 
rate is the rate at which new protons are introduced into 
t'e rotating planetary raagnetosphere (assuming it exists 
for Saturn) . if charge exchange is the dominant lifetime 
process for hydrogen , the introduction of a new proton 
requires the loss of an old piroton, so that only the energy 
distribution, not the number of protons is changed. A net 
production of protons will however result from elecron impact 
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RADIAL DISPLACEMENT FROM JUPITER (Rj) 


Figure 29 

lite Proton Source. Model results for the radial distributions of 
OSS rate of hyarogen atonis frciu the extended atmospheres of Ztoialth 
ede and Callisto by charge exchange processes are shown. The tota 
ation of each cloud and its overall loss rate are indicated. See 
ext for discussion. 
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ionization of the extended hydrogen atmospheres, with a rate 

-2 

of approximately 10 times the charge exchange rate. The 
neutral atoms released from the charge exchanged process are 
on gravitational escape orbits from the planet if initially 
located outside of a critical radius. Those neutral hydrogen 
atoms initially located inside this critical radius are on 
gravitational bound orbits of the planet and upon re-ioniza- 
tion may enhance the proton density locally. For Jupiter and 
Saturn the critical radius (located where the magnetospheric 
rotational and orbital escape velocities are equal) is 2.83 
and 1.86 planetary radii respectively. Such an enhancement 
is then only possible for the hydrogen cloud of Amalthea 
where the satellite orbit is 2.54 Jupiter radii. 

The impact of the satellite protons on the Jovian mag- 
netospheres needs to bo assessed by compairing this source 
with those provided by the solar wind and the planet. In 
the case of Titan and Amalthea, the net production of protons 
could be largo enough to inflate the magnetosphere locally 
and create unique planetary plasraaspheres . 

4.3.2 Sodium Ion Source 


Models for the sodium cloud of lo (Smyth and McElroy, 
1978) suggest that neutral sodium is emitted from the inner 
hemisphere with a mean initial velocity of 2.6 km/sec and 
an average hemispherical flux of 10 atoms cm~ sec” . 

Atoms are lost from the extended sodium atmosphere by 
electron impact ionization with an effective lifetime of 
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between 15 and 20 hours, -rhe radial distribution for the 
loss of sodium atoms from the cloud for a 20 hour lifetime 
is shown in Figure 31. The loss rate has a maximum value 
of 2.6 X 10^^ atoms per 0.1 Jupiter radii per second at the 
satellite position and the cloud is localized between about 
4.2 and 8 Jupiter radii. 

For a steady state cloud, Figure 31 gives the radial 
distribution for production of sodium ions in the Jovian 
magnetosphere Once formed, the heavy ions are swept away 
by the rapidly rotating planetary field and become part of 
the local plasmasphere. Radial diffusion of these heavy 
ions will then follow and way distribute the ions in such a 
way so as to produce centrifugal distortion of the Jovian 
magnetic field. 

4.4 Modeling lo's Sodium Line Profile Data 

A model capable of calculating the sodxum Dl and D2 
line prof iles generated by lo' s extended sodium atmosphere 

when observed through a given viewing aperature has been 
developed. The positions and velocities of the satellite, 
sun and Jupiter may be selected to match a given observation 
The modeling parameters are the absolute flux, the distribu- 
tion of flux on the emitting satellite surface or exosphere, 

the initial emission velocity distribution and the lifetirae 
of sodium atoms in the Jupiter environment. Preliminary 
efforts have been directed to understanding the effect that 
these modeling parameters have upon the line profile shape 
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RADIAL DISPLACEMENT FROM JUPITER (Rj) 


- — Figure 31 

lo. Model results for the radial 
drstribution of the loss rate of sodium atoms from lo’s extended 
atmosphere by electron impact ionization processes are shown. 

The total atom populatxon and overall loss rate of the cloud are 
indicated. See the text for discussion. 



for a variable satellite phase angle. 

Comparison of modeled line profiles with existing obser 
vations is being used to limit and define the range of the 
modeling parameters. In addition, a cooperative effort has 
been established with observers who are making specific 
measurements needed to test for the presence or absence of 
particular line profile features generated by model calcu- 
lations. Progress to date has been most encouraging, but a 
considerable effort will be required before a consistent 
and accurate interpretation of the sodium line profile data 
is available. Efforts to accomplish this line profile 
modeling goal will be pursued in continuing research work 
funded by NASA . 



CHAPTER V 


CONCLUDING REMARKS 


Exploratory and Prototype modeling of hypothetical 
extended hydrogen atmospheres of the outer satellites 
Amalthea, Ganymede, Callisto and Titan has been presented. 
These results, given in Chapter IV, were calculated using 
our cloud model discussed in Chapter III, and a variety of 
atom lifetime and emission conditions were considered. For 
the most probable conditions, discussed in Chapters II and 
III, the extended hydrogen atmosphere of Amalthea forms a 
tightly-bound partial toroidal-shaped cloud about its 
planet, whereas Ganymede, Callisto and Titan - having much 
U^s^ter lifetimes — form rather large, complete and nearly 
symmetric toroidal-shapod clouds. Model results suggest 
that the partial toroidal hydrogen cloud for Amalthea, 
tentatively detected by the UV instrument of Pioneer 10 
(Judge et al . , 1976), could be explained if the satellite 
atom escape flux is of order 10^^ cm ^sec~^. Results also 
suggest that the Voyager spacecrafts should be able to detect 
the Lyman-a emission from the extended hydrogen clouds of 
Ganymede, Callisto and Titan, assuming most probable condi- 
tions, and for Ganymede and Titan even under less favorable 
conditions . 

The initial ion sources, created when cloud atoms of the 
extended hydrogen atmospheres of Amalthea, Ganymede, Callisto 
and Titan are ionized, were calculated. These sources, 
together with the heavy ion source created by the sodium 
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cloud of lo were presented in Chapter IV. The impact of 
these ions upon the planetary magnetospheres of Jupiter and 
Saturn (assuming it exists and includes Titan's orbit) may 
be significant and warrants detailed investigation. 

Progress to date in modeling lo's sodium line profile 
data is discussed in the last section of Chapter IV. Using 
our fully-developed cloud model, initial efforts have been 
focused upon limiting and defining the range of modeling 
parameters. A variety of sodium observations will be used 
to accomplish this goal to be pursued in our continuing 
research work funded by NASA. 

Results presented in Chapter IV should be particularly 
useful for preliminary interpretation of data acquired in 
1979 by the Voyager I and Voyager II spacecrafts encountering 
Jupiter njid the Pioneer .11 .spneocraft encountering Saturn. 

In addition to these particular model results, information 
presented in Chapter III may be used to estimate the spatial 
extent of any extended satellite atmosphere and to estimate 
the atom emission and lifetime conditions from spatial extent 
data measured by spacecrafts. Analysis of such spacecraft 
data should yield important clues as to the nature- of the 
local or bound satellite atmospheres, the character of the 
electromagnetic interaction of satellites with their plane- 
tary magnetospheres and the concentration, energy density 
and spatial distribution of charged particles in the mag- 
netospheres. 
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